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Abstract. We introduce an effective model for e g electrons to describe three- 
dimensional perovskite (Lai-zSr^MnOs and Lai-^CazMnOs) manganites and study 
the magnetic and orbital order on 4 x 4 x 4 cluster using correlated wave functions. 
The model includes the kinetic energy, and on-site Coulomb interactions for e g elec- 
trons, antiferromagnetic superexchange interaction between 5 = 3/2 core spins, and 
the coupling between e g electrons and Jahn- Teller modes. The model reproduces the 
experimentally observed magnetic order: (i) A-type antiferromagnetic phase in the 
undoped insulator LaMnO*3, with alternating e g orbitals and with small Jahn- Teller 
distortions, changing to a conducting phase at 32 GPa pressure, and (ii) ferromagnetic 
order in one-eight doped Lay/gSr^gMnOs and in quarter doped ha 3 /^Sv 1 / 4 M-'a03 com- 
pounds. For half-doped La^Ca^MnOa one finds a competition between a ferromag- 
netic conductor and the CE insulating phase; the latter is stabilized by the Jahn- Teller 
coupling being twice larger than for the strontium-doped compound. Altogether, there 
is a subtle balance between all Hamiltonian parameters and the phase diagram is quite 
sensitive to the precise values they take. 
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1. Introduction 

Much effort was put into understanding what is going on in doped perovskite manganise 
oxides both on theoretical [IHT3] and experimental P3H26] side. Early expectations 
that simple physical mechanisms, such as involving for example only Jahn-Teller 
(JT) interactions and Hund's exchange but without Coulomb interaction, or pure 
electronic Coulomb interactions while neglecting Hund's exchange, could explain the 
phase diagram and the phase transitions observed situations in the manganites were 
not confirmed [HE] • Also relatively early it was realized that the core t2 9 electrons, even 
if dynamically passive, play an essential role for the stability of magnetic and orbital 
order [HE]. Altogether, the evidence has accumulated that there is a delicate balance 
between several competing physical mechanisms which lead to rather complex phase 
diagram of the perovskite manganites [HIUKIS] . 

Investigation of doped manganites is challenging as the Jahn-Teller distortions 
may lead to the charge order which will also favour particular orbital order [T|l2~T]. 
To make the complex situation in doped manganites tractable in the theory, several 
theoretical papers were focused in past on single-layer and double-layer manganites 
as the description of quasi two-dimensional (2D) systems was simpler than the one 
necessary for three-dimensional (3D) doped perovskite manganites (see for example [11] 
and the references therein). Realistic 3D systems were more difficult to tackle at least on 
a better basis than simple mean-field type approaches. Due to strong local correlations, 
however, such approaches were not quite appropriate to provide a better insight into 
the competition between different possible physical mechanisms present in these strongly 
correlated electron systems. The present paper is one attempt (among the others) to 
fill up this gap. 

Similar to model studies [5|fTT] . it has been also recognized in the electronic structure 
calculations performed within the local density approximation (LDA) that the local 
interactions are important and one has to use either LDA+?7 or LDA+DMFT |12,28j, 
i.e., take into account also the dynamical aspects in the dynamical mean field theory. 
These calculations confirmed the earlier point of view [5JI29] that electron-electron 
interactions and the Jahn-Teller terms support each other and both are necessary to 
reproduce the observed orbital order in the undoped compound. In particular, by the 
studying the effect of pressure it was concluded that LaMnOs is not a Mott-Hubbard 
insulator, but rather a Jahn-Teller insulator, however, it has been recognized that the 
correlations play also an important role [28] . 

The paper is organised as follows. First, we introduce (in section [2]) a realistic 
model for e g electrons in the perovskite manganites which includes the electron Coulomb 
interactions, the coupling between e g electrons with s = 1/2 spins and t 2g core 
5 = 3/2 spins, superexchange interaction, and the interactions of e g electrons with 
local distortions due to the Jahn-Teller terms (section 12.11) . In this section we also 
present the method to determine approximate ground states using a 4 x 4 x 4 cluster. 
We treat electron correlation effects in the ground state beyond the Hartree-Fock 
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(HF) approximation using the so-called local approach (section 12.21) . The numerical 
results are presented and analysed in section [3j Here we address a question which 
values of the parameters of the effective Hamiltonian are appropriate to account for 
the experimentally observed situation and present results obtained for the undoped 
LaMnOs, and for Lai_ x Sr x Mn03 at doping levels x = 0.125, 0.25 and 0.50. Finally, 
section H] contains a short summary of the paper and general conclusions. 

2. The model and calculation method 

2.1. The model for e g electrons 

We study strongly correlated electrons in undoped and doped 3D perovskite manganites 
Lai_3.Sr3.MnO3 (or Lai-^CajMnOs) using effective model describing only Mn sites 
renormalised by surrounding oxygens. At each site we consider Wannier orbitals of 
the e g character composed out of manganese and oxygen orbitals [12J. One of them 
is filled at Mn 3+ ions, and both are empty at Mn 4+ ions in a doped system. In both 
cases the t% g orbitals on Mn ions are occupied by three "core" electrons with total spin 
5 = 3/2 (treated as frozen and classical). Thus the active electrons here are only e g 
electrons. We investigate the system by using a model Hamiltonian 

H = H kin + H int + H spin + iJ JT (1) 

which consists of the kinetic energy, the on-site Coulomb interactions, spin interactions, 
and the Jahn-Teller term. This model and/or the essential parts of it, were studied 
earlier by several groups, see refs. [T|l?tl8l [TTHT3l [T9 |l3lM?3"] . 

The kinetic part Hkin is expressed using a local e g orbital basis at each site, which 
reads in short notation, 

\z) = (3z 2 - r 2 )/V6, \x) ee (x 2 - y 2 )/V2. (2) 
Using this local e g orbital basis at each site one finds anisotropic phase- dependent 
hopping 0[EEJ[T2] 

H iin — — -t ^2 { (3d\ X(T dj XfT + d\ m dj Za ) ± V3(dl xa dj Za (3) 

{ij}\\ab,cr 

{v}\\c,o- 

Here d\ are creation operators for an electron in orbital fj, — x,z with spin a =f, \, at 
site i. The {i, j} runs over pairs of nearest neighbours and gives two contributions for 
each bond (i, j) ; ± is interpreted as plus sign for the bond (i,j) being parallel to the 
crystal axis a and minus for the bond (i, j) parallel to the axis b. (The crystallographic 
axis c is assumed to be aligned with Cartesian z-axis.) 

The i?i nt and H spin terms follow from the Coulomb interactions within a degenerate 
d band [31]; these terms are: 

5 

Hint = Up + (^0 - 2 Jh ) ^2 n ix n iz, (4) 

ifi i 
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(7) 



Here H mt describes the charge interactions within the e g subsystem and plays a role 
in doped manganites. The spin interactions H spin include the leading part of Hund's 
exchange Jh, both between two e g electrons at the same site, and between s = 1/2 spin 
of a single e g electron and the 5 = 3/2 core t 2g spin at each site. On-site intraorbital 
Coulomb interaction is denoted as Uq, and the interorbital Coulomb interaction is a 
linear combination of this term and Hund's exchange Jh [31] ■ Note that we include 
here only the leading Ising part of Hund's interaction terms, which is an approximation. 
However, within the present method using the correlated wave functions which are 
obtained by correcting the Hartree-Fock (HF) approximation, the remaining transverse 
terms do not contribute, compare refs. ) - The term oc J' is the antiferromagnetic 

(AF) superexchange generated by charge excitations of t 2g electrons, which lead to AF 
Heisenberg interaction between frozen core t 2g electrons |8|H2]I29]. This frozen-core 
approximation for S = 3/2 core spins works well for the description of the ground state 
properties as shown by earlier studies, see e.g. refs. [TT|l35]. 
Finally, the Jahn-Teller Hjt part is 



#JT = Yl {^ JT (^foiz + + Q^i + Q^i) (8) 

i 

+^Wu + Ql i + Ql i )\- (9) 



2 

and includes three different Jahn-Teller modes {Qu, Qn, Q31} at each site. The operators 
Tf in ifjx are constructed for r = 1/2 orbital pseudospin for e g electrons at site 

T i ~ ^ ^ {djxo-diza + d iza di xo ) , 

(J 

T i = ^2(4xAxa ~ 4zAzv), (10) 

(7 

while Qu, Qu and Q^i denote the active Jahn-Teller deformation modes of the 
i— th octahedron. (For simplicity the harmonic constant of the isotropic Jahn-Teller 
(breathing) mode Q\ is assumed to be double with respect to those corresponding to 
Qi and Q3 unsymmetric modes as discussed in ref. pQ). The breathing mode is quite 
often neglected in the effective models for the perovskite manganites, and plays no role 
in the undoped system. As we are dealing here with doped systems, we decided to keep 
it to investigate its consequences. 

The crystal field splitting between z and x orbitals (j2j) was assumed to be zero 
(this approximation seems reasonable in quasi-cubic 3D manganites [13] though some 
authors [12] assume instead small finite values). Below we describe the numerical 
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simulatons performed on finite 3D clusters with the above model ([T]) The variants of the 
described model reproduced the observed sequence of magnetic phases for increasing 
hole doping in 2D monolayer and bilayer manganites [32,33]. 

2.2. Cluster approach including electron correlations 

We studied 4x4x4 clusters with periodic boundary conditions (PBC) filled be different 
number of e g electrons, corresponding to the undoped LaMnOs system, and to systems 
with doping one-eight (x = 0.125), quarter-doped system (x = 0.25) and half-doped 
system (x = 0.5). All calculations were performed at zero temperature (T = K). 
First, the calculations within the single-determinant HF approximation were performed 
to determine the ground state wave function |$o)- Also the energetic distance A between 
the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular 
orbital (LUMO) (HOMO-LUMO gap) was extracted at this step. 

In the next step each HF wave functions was independently modified to improve 
the energy and to include the electron correlations using the so-called local ansatz (for 
details see refs. [30-33J). We used the ansatz for the correlated ground state, 

|M/)=exp(^r /j O J )|$o), (11) 

i 

where {rji} are variational parameters, and {Oj} is a set of local correlation operators. 
These operators include the leading local density-density correlations in the present 
model for the perovskite manganites, in analogy to the 2D layered manganites, see refs. 
[32j[33] . These local operators used in the present model correspond to the subselection 
of the most important two electron excitations within the ab initio configuration- 
interaction method. The variational parameters {rji} were found by minimising the 
total energy 

Etot -^W (12) 

In this way the correlation energy, 

-Ecorr = E tot — -EfjF, (13) 

was obtained. 

Coming to technical details, in the first step HF computations were performed for 
each considered electron filling, starting from one of several different initial conditions 
(about a thousand for each set of Hamiltonian parameters), i.e., from predefined 
(some symmetric but mostly random) charge distribution, spin configuration, the latter 
selected from eight predefined patterns of core t 2g spins (see figure 1 of ref. [33]), and 
from several predefined sets of classical variables {Qu, Q 2 i, Qzi}- F° r eac h fixed set of 
the starting parameters and starting initial conditions on convergence we obtained a 
HF wave function which was a candidate for the ground state wave function. This self- 
consistent procedure was performed to provide energy minimum also with respect to 
the {Qu, Q21, Qu} classical variables 
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Figure 1. A simplified phase diagram in (U, Jh) plane for undoped (x — 0) LaMnC>3. 
G refers to an ordinary Nccl antifcrromagnct (G-AF phase), A to A-AF phase with 
FM planes and AF coupling along the direction perpendicular to them, and F stands 
for the FM phase. Parameters: t = 0.4 eV, g JT = 2.2 eV A -1 , K = 13 eV A" 2 and 
J' = 3 meV. 



After completing the HF computations we performed correlation computations (as 
a second step for each investigated wave function) obtaining the total energy ( |T2l) and 
the correlation energy ( !T3|) (for more details see refs. [30-32J). After obtaining the total 
energy for one configuration, we repeated the entire procedure from the beginning, i.e., 
we took another (second) set of HF initial conditions and repeat all computations to 
obtain the second candidate for the ground state. Similarly for the third, fourth, • • •, 
etcetera, set of initial conditions. Finally, the resulting set of total energies was inspected 
and the lowest one was identified as a good candidate for the true ground state. 

3. Numerical results 

As the parameters of the model (CQ) are known only with certain accuracy, the 
computations were repeated for many sets of the Hamiltonian parameters. The values 
of the parameters which one should consider as realistic for La 1 _ 2 .(Sr,Ca) x Mn03 were 
extensively discussed in the literature (for a review in context of the present paper 
see refs. [HEED [32l> [33]). Following the standard parameter sets [32], [33], we assumed 
that to — 0.4 eV and the intraorbital Coulomb element is limited in the range 
5 < U/t < 12. The values of Hund's exchange are larger than t [9l [29], [38] , and 
we considered to < Jh < 2to- The Jahn- Teller constant K was fixed as K = 13 eV 
A~ 2 [U[2][30][32l[39]. Not much is known about the value of the coupling constant gjT, 
therefore following discussion in ref. p] we assumed that 2 eV A~ x < < 3.7 eV A -1 . 
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Figure 2. Charge distribution, magnetic, and orbital order of e g electrons in the A- 
AF phase of LaMn.03, as obtained in undoped 4x4x4 cluster with PBC. charge and 
alternating x, z orbital order for e g electrons Only single ab planes are shown. Core ti g 
spins are not shown. Left side panels — At each site the circle diameter corresponds 
to total e g on-site charge; the arrow length to the e g spin; the horizontal bar length to 
charge density difference between x and z orbitals (longest bar to the right corresponds 
to pure x, to the left to pure z, zero length to half by half combination). All these 
values are expressed in approximate proportionality (as generated by graphic package 
of latex) to nearest neighbour site-site distance which is assumed to be unity. Upper 
panels show situation on even numbered ab planes, lower panels correspond to odd 
numbered ab planes. Right side panels — Jahn-Teller distortions for the same even 
and odd ab planes: Qn (bars drown slightly to the left of each site) and Q%i (bars to 
the right of each site), and isotropic Q\ (breathing) mode in between. For more clarity 
Qi and Q3 bars are enlarged by factor of 2 compared to Q\. In the present case all 
Qu take the value of -0.08 A, the largest Q^i is equal to 0.10 A, the largest Q-a is 0.14 

A. 

And, finally, following recommendation of refs. [13|l29] we use the experimental value of 
the Neel temperature in CaMnOs [ID], where e g electrons are absent and the isotropic 
superexchange between core t 2g spins follows solely from the excitations of t 2g electrons, 
to fix J' = 3 meV. Thus we have three fixed and three free parameters. 

The aim was to identify the set (or sets) of these three free parameters which 
lead to good agreement of the model predictions with experimental data. This task 
turned out to be rather difficult especially with respect to Jahn-Teller parameter g JT . 
Actually almost all values of gj^ either yield conducting state for both x = 0.125 and 
x = 0.25 doping levels, or they yield an insulator, also for both above doping levels. This 
contradicts the experimental situation as Lai-^Sr^MnOa manganite is weakly insulating 
for x = 0.125 and conducting for x = 0.250 [H]. This insulator-to- metal transition 
turned out to be very useful as it allowed us to fix all three free parameters. The values 
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Figure 3. HOMO-LUMO gap A for increasing pressure p in the undoped LaMnOs: 
results of the calculations are shown by circles, solid line is a guide to the eye. The linear 
dependence of to(p) and K(p) on pressure is assumed. Other Hamiltonian parameters 
take standard values, see the caption of figure 2. 

of the parameters which can bring out the agreement with the experimental data form 
a small isolated region within the parameter space, and finding them required making 
hundreds of scans through the parameter space. We have found that only for U/to — 8, 
Jh = 0.7 eV and #jt = 2.2 eV A -1 we are able to reproduce the experimental results. 

3.1. Undoped LaMn0 3 

The insulating phase in LaMnOs was reproduced by the present calculations. For the 
magnetic order we obtained rubust A-AF structure, which has ferromagnetic (FM) order 
in ab planes, coupled by AF interaction in the c direction. This phase is almost generic 
and can be reproduced on the phase diagram for many sets of Hamiltonian parameters. 
Morover, the correlation energy is small (see Table 1) and does not change the magnetic 
and orbital order obtained from the HF wave functions. To give an example, let us 
show one 2D section of our simplified phase diagram. Note that as long as we ignore 
the global structural distortion of the lattice below that structural phase transition, all 
cubic directions are equivalent and A-AF order shown in figure 1 applies also e.g. to 
FM order in ac planes accompanied by the AF coupling along the b direction. 

As already mentioned the experimental data on conductivity fix the electron 
interaction parameters U/t Q = 8 and Jh = 0.7 eV, as well as #jt = 2.2 eV A -1 . For these 
parameters the ground state ordering for LaMnOs is shown in figure 2. The Jahn- Teller 
Q2 and Q3 distortions measured in ref. [23J (at T = 200 K) are only roughly in agreement 
with our results. (Note that for simplicity the cubic Jahn- Teller terms [71 [13] in ifjx 
were omitted. Therefore, the Jahn- Teller distortions we obtained are only estimations 
rather than the precise values). 

We have found a large HOMO-LUMO gap A for LaMnOs which is an indicator 
that the obtained A-AF phase is insulating, see figure 3. It is known that at pressure of 
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Figure 4. HOMO-LUMO gap A for increasing Jahn- Teller coupling <7jt- Circles, 
squares and crosses correspond to x — 0, x = 0.125 and x = 0.25, respectively. 
Parameters: t = 0.4 eV, U/t = 8, J H = 0.7 eV, J' = 3 meV. 

32 GPa LaMnOs becomes conducting [42J. In order to investigate the effect of pressure 
we followed heuristic argumentation presented in ref. [43] - it has been argued there 
(see figure 2 in ref. [43]) that the effect of 32 GPa pressure can be reproduced in a 
model similar to our effective model (JTJ) when to and K Hamiltonian parameters are 
renormalised in the following way: to — > 1.4t and K — > 3.3K. Following ref. [13], the 
other parameters (of the model) were left unchanged. We note that this assumption 
is somewhat unrealistic as at least the superexchange J' would also increase under 
pressure, but this effect is not important as long as we investigate the ground state 
and the magnetic structure does not change. Uo and other effective parameters of the 
Hamiltonian should also show some limited variation under increasing pressure. 

Inserting the new renormalised values of hopping and Jahn- Teller constant into 
our computations we found that the high pressure phase is indeed a conductor and we 
could identify the conducting phase as perfect charge-homogeneous, FM phase with zero 
Jahn- Teller distortions and with perfect half-to-half (x 2 — y 2 ) and (3z 2 — r 2 ) orbital order 
on each individual site. This orbital liquid state occurs in the doped FM manganites 
under normal pressure [44] . 

To provide a better evidence to such a prediction we plotted HOMO-LUMO gap 
A for increasing pressure. For simplicity it was assumed that to{p) and K(p), where p 
denotes pressure, vary linearly from i an d K at p = to 1.4t and 3.3K at p = 32 
GPa. Anyway, the transition from insulating to conducting phase can be deduced from 
figure 3 and it occurs around p ~ 30 GPa where the gap becomes small enough, A ~ 0.1 
eV. Note that in the present finite cluster one finds always a gap between HOMO and 
LUMO, but we have estimated that A ~ 0.1 eV corresponds to a metallic ground 
state in the bulk. To close the discussion on LaMn0 3 under pressure let us mention 
ref. [45] where magnetic and orbital order under uniaxial strain (only JT modes were 
pressure-renormalised) was studied by using similar (to our) Hamiltonian. 
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Table 1. Hartree-Fock energies E to t and total ground state energies -Erf (fl2|) - both 
in eV, as obtained for 4 x 4 x 4 cluster versus the doping x. Parameters of the 
Hamiltonian: t = 0.4 eV, U = 8t , J H = 0.7 eV, K = 13 eV/A 2 , g = 2.2 eV/A, 
J' = 3.0 meV. 

X EhF ^tot 

0.000 -96.016 -96.114 

0.125 -88.217 -88.315 

0.250 -79.374 -79.521 

0.500 -58.666 -58.896 



3.2. One- eight- doped Lao^Sro. 125 MnOs 

In one-eight doped manganite one finds an insulating phase with FM order, both in the 
HF and for the correlated wave functions. The correlation energy is small, see Table 
1, and does not change the magnetic and charge order in the system. The breathing 
mode Jahn-Teller distortions are small and uniform, Qu f» —0.08 A. This phase is 
characterized by a rather interesting charge and orbital distribution, with a distinct 
alternation in charge density between even and odd planes. Consequently, this phase is 
an insulator, in agreement with the experiment [41] . though a weak one as indicated by 
rather small HOMO-LUMO gap shown in figure 4 (curve in the middle, squares). Note 
that gap is reduced by a factor close to 3 from the value obtained for LaMn0 3 . 

The alternation of e g electron density occurs between even and odd planes. First, in 
even numbered planes there are somewhat less electrons and the orbitals are of (x 2 — y 2 ) 
character, Q21 ~ and Qsi take the value of —0.1 A. These distortions favor the x 
orbitals. However, there are also some exceptions, i.e., occasionally on a site we find the 
(3z 2 — r 2 ) orbital occupied and Q 3i ~ +0.1 A. The charges are distributed in a slightly 
non- homogeneous way, with the site to site variations of charge being at most 0.2 e. 
Second, in the odd numbered planes one finds a little bigger average charge density per 
site than in the even planes and almost homogeneous charge distribution. In contrast 
to the even planes, one finds predominantly (3z 2 — r 2 ) orbitals occupied. The Jahn- 
Teller distortions are Q 3i « ±0.1 A and Q 2i w ±0.1 A. A gain, we have seen also some 
exceptions, with x orbitals occupied at a few sites and the corresponding modified Jahn- 
Teller distortions, similar to those in the even planes. In fact, this result shows that 
sites with differently occupied orbitals can proliferate between the even and odd planes 
and form charge and orbital defects in each plane. This shows that the inhomogeneous 
charge distribution is a generic feature in this range of doping. 

3.3. Quarter-doped Laoj^Sro^MnOz 

Also for quarter-doped manganites, such as Lao.75Sr a 25Mn0 3 , one finds FM order and 
weakly non- homogeneous charge distribution with the present parameters (see figure 5). 
These results are similar to the ground state found for x = 0.125 — we have found again 
nonequivalent alternating planes, as described below. The correlations do not change 
this result, and are again weak (Table 1) as the local moments have formed and the 
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Figure 5. Magnetic A-AF, charge and orbital electron density for e g electrons as 
obtained in quarter-doped manganites, x = 0.25. Four single ab planes (perpendicular 
to the c axis) extracted from 4x4x4 cluster are shown. Core tig spins are not shown; 
their directions agree with the directions of e g spins. The meaning od symbols and 
parameters are the same as in figure 2. 

main effects are already captured in the HF approach. 

On the one hand, within even numbered planes in the ground state one has weak 
charge order in form of a stripe phase, with vertical or horizontal charge-minority lines 
alternating with charge-majority lines (charge difference of about 0.2e). Charge minority 
sites have mixed half to half x and z orbital occupancy and negligible Jahn-Teller 
distortions, which is characteristic for the orbital liquid in the metallic state [44]. The 
situation is quite different along the charge majority lines — one finds there half of these 
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sites with x orbitals and with negligible {Qu,Q2i} distortions, and with Q 3i —0.13 
A. The other half of these sites carries the mixture of z and x roughly in proportion 
2:1 and negligible {Qu,Qsi} distortions, whereas Q21 ~ —0.1 A. These sites with higher 
electron density (i.e., those with x orbitals and those with 2:1 orbitals z and x) alternate 
along the charge- majority lines. 

On the other hand, on odd numbered planes the orbitals and Jahn-Teller distortions 
are much the same as the ones on even numbered planes (as described above). However 
there is a difference — along charge- minority lines (each second line in the cluster) we 
did not find homogeneous charge density but, instead alternating charge- minority and 
charge- majority sites. 

The FM phase at x — 0.25 doping should be a metal which follows from quite 
small HOMO-LUMO gap (A ~ 0.1 eV), see figure 4. This result agrees with the phase 
diagram for the perovskite manganites [l,46j. 

3.4- Half-doped LaQ^Sr 5 MnOs — stability of the CE phase 

Surprisingly, for half-doped Lao.sSro.sMnOs one finds again almost the same ground 
state as for x = 0.25, namely FM spin order, accompanied by weakly non-homogeneous 
charge distribution and a conducting phase (HOMO-LUMO gap is about 0.12 eV). The 
ratio of x to z orbitals on each site is half to half, while the Jahn-Teller distortions 
are very small. These findings reproduce the experimental situation in Lao.sSro.sMnOs, 
where no CE phase was observed. 

On the contrary, the CE phase was found in half-doped Lao.sCao.sMnOs [17JHH]. 
A somewhat unprecise but still common point of view is that a manganite doped with 
strontium is more metallic than the manganite doped with calcium (more insulating) 
which in turn might be attributed to weaker coupling to the lattice by the Jahn-Teller 
interactions in strontium compounds and a bigger role of these interactions in calcium 
compounds. Indeed, the FM phase was found insulating in a broader range of doping in 
Lai-ajCa^MnOa than in Lai^Sr^MnOs [19] . If this is the case then our effective model 
should provide a different prediction for the ground state when g JT is larger. 

When we discussed the precise values of Hamiltonian parameters following the 
discussion presented in ref. [1] we assumed 2 eV A" 1 < gjx < 3.7 eV A -1 . Then for 
strontium manganite Lai_. r Sr.j,Mn03 we found that gjx = 2.2 eV A -1 . Now for the sake 
of the ongoing discussion for calcium doped manganite we will accept gjT = 3.7 eV A -1 , 
the biggest value which is allowed. This change of the unique parameter immediately 
provides us with the expected result. Now the spin order in the ground state is indeed in 
form of FM zig-zags with AF order between them, characteristic for the CE phase, and 
the charge order is of checkerboard type (figure 6). The CE phase is quite robust, it is 
identified at the HF level and the correlations when added do not change this result. This 
is in contrast to what happens in bilayer half-doped manganites where the correlations 
stabilize the CE phase [33]. We note, however, that the magnetic moments are larger 
in the 3D system, and therefore the role of electron correlations in the magnetic states 
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Figure 6. Magnetic CE, the ordered charges and fractions of occupied x and z orbitals 
for e g electrons in half-doped Lao.5Cao.5Mn.O3. The standard Hamiltonian parameters 
(which apply to strontium-doped substance) were assumed as in the caption of figure 
2 with one exception, namely we adopted a larger Jahn-Teller coupling gjx = 3.7 eV 
A -1 . Two panels correspond to two neighboring nonequivalent ab planes — the order 
found in these planes is repeated along the c axis in the planes number 3 and 4, not 
shown. Only e g spins are shown (T2 g spins are parallel to them at all sites) . Note that 
the magnetic CE zig-zag-like order refers to the total spin thus it is not properly visible 
on the figure where only e g spins are plotted. The legend is the same as in figure 2. 

such as the CE phase is reduced. 

3. 5. Switching off John- Teller distortions 

The results presented up to now indicate the importance of Jahn-Teller interactions. 
Due to them local distortions are created at each doping level, and they influence both 
the electronic structure and the magnetic properties. It seems that the value of gjx 
coupling is the primary factor which determines whether the 3D perovskite system is 
a conductor or an insulator. Therefore it seems reasonable to perform additional test 
computations switching off the Jahn-Teller interactions. 

We performed such computation for Lai.^Sr^MnOs where x = 0,0.125,0.250 and 
gjT = 0. The resulting ground states were all the same: FM, charge homogeneous orbital 
liquid phase jB], with half to half ratio of x to z orbitals at each site and conducting 
(HOMO-LUMO gaps we found were all about 0.1 eV). This confirms the above point 
of view that the Jahn-Teller distortions play a crucial role for the charge, magnetic and 
orbital order in doped perovskite manganites. 
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4. General conclusions and final remarks 

We considered the effective Hamiltonian for e g electrons which takes into account 
their kinetic energy with anisotropic phase dependent hopping oc t, on-site Coulomb 
interaction U, Hund's exchange coupling Jh, their coupling to frozen core t2 9 electrons 
oc Jh, the Heisenberg AF superexchange oc J' between core T 2g spins, and the Jahn- 
Teller interactions with lattice distortions. Due to the complexity of the problem, we 
could not provide a precise proof that this model is sufficient for the physical properties 
of the perovskite manganites, but rather we report on results of computations which 
demonstrate that all these interactions enter on equal footing and seem to be essential 
for correct interpretation of the experimental data. Indeed, we have performed also 
calculations with several sets ignoring one or the other interaction term, and the results 
were very unsatisfactory in all cases, i.e., such simplified models were not able to 
reproduce the experimental observations. 

We would like to emphasize that the realistic Jahn-Teller coupling plays a rather 
special and a very important role in the manganites as it controls the conductivity of 
the perovskite 3D systems. We have presented evidence that this interaction generates 
local distortions in the entire investigated range of hole doping 0.125 < x < 0.50 which 
lead to charge inhomogeneities, defects in the orbital occupancy, and influences the 
magnetic order. Such states have been seen before in the experimental data on doped 
manganites [50] and were also suggested to follow from the theoretical models |51j . 
Unfortunately, reliable studies of systems with charge disorder and local distortions 
are very difficult to perform, but the present results demonstrate that these states 
might decide about the colossal magnetoresistance and select one or the other type of 
magnetic order. Particularly at half-doping (x = 0.5) we have found a subtle balance 
between the FM phase and the CE phase, the latter stabilized by large Jahn-Teller 
interactions. This is not surprising as the nearest neighbor JT coupling might not be 
sufficient to stabilize the observed zigzag FM chains in the CE phase for the realistic 
parameters. The mechanism invoked in |52j to stabilize the CE phase is subtle and 
employes the cooperative JT interaction between next-nearest Mn 3+ neighbors mediated 
by the breathing mode distortion of Mn 4+ octahedra and displacements of Mn 4+ ions. 

As found before for monolayer and bilayer doped manganites, see refs. [301132] . we 
have obtained many metastable ground states in the 3D system which are possible in 
principle. Such states are not discussed in detail in the present paper, but we just 
point out that they represent local energy minima and could be also realized at finite 
temperature. Therefore, the search for the best candidate for the true ground state 
requires examination of numerous phases with different kinds of nonhomogeneous charge 
distribution. Therefore, we suggest as a general conclusion that nonhomogeneous charge 
distribution is a generic feature of the doped perovskite manganites. It follows that the 
theoretical predictions which are based on models assuming just a few a priori selected 
candidates for the ground state (especially the ones with high symmetry of charge or 
orbital order) can not be fully conclusive unless they are based on reliable experimental 
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data. 
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